function [dftoYVec,deficitVec,sampVec] = deficitSequence( Nquarters , dftoYInitial , rho)
sampVec=1998:0.25:2008.00;
sampVec=sampVec(1:Nquarters); 
deficitVec=fliplr( linspace(-0.0589,-0.0164,Nquarters) + 0.0142 ); 
dftoYVec=zeros(Nquarters,1); 
dftoYVec(1)=dftoYInitial - deficitVec(1)/(1-rho); 
for ii=2:Nquarters; 
   dftoYVec(ii)= dftoYVec(ii-1) - ( deficitVec(ii)-deficitVec(ii-1) )/(1-rho)...
       -deficitVec(ii-1); 
end 
